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Chapitre 1 

Introduction et définitions 

1.1 Définitions 

1.1.1 Prologue 

Les mathématiques ont été et sont appliquées, utilisées pour contrôler les relations des hommes 
avec son environnement. Parmi les travaux qui ont marqué l'histoire, nous citons: 

■ Archimede (-285, -215 AJC): loi de la statique. 

• J. Kepler (1571, 1630): mouvement des planètes et sections coniques. 

• I. Newton (1643, 1727): loi de la dynamique et satellites artificiels. 

• L. Euler (1707, 1783): travail sur les mâts des bateaux. 

• Gauss (1777, 1852): méthodes des moindres carrés. 

• B.Riemann (1826, 1866): physique mathématique. 

• H. Poincarré (1854, 1912): mécanique celestre. 

Depuis la deuxième guerre mondiale, les applications des mathématiques s'étendent à tous les 
secteurs d'activité et ceci grâce aux ordinateurs et leur rapidité de calcul: on peut citer les prévisions 
météologiques, les techniques d'exploration de la terre en géophysique comme la sismique, les 
techniques de recherche opérationnelle en gestion, etc.. 

Un cours de "Mathématiques Appliquées" qui exige de faire des mathématiques avec un ordi- 
nateur, se trouve alors justifié. Si les autres disciplines de mathématiques comme (Algèbre linéare, 
Analyse vectorielle... ) sont connues depuis très longtemps, le terme analyse numérique n'a com- 
mencé à être utilisé qu'après la fondation en 1947 de "Institute of Numerical Analysis" à l'Université 
de California de Los Angelos. 

1.1.2 Définition de l'Analyse numérique (A.N). 

Il sera difficile de donner une définition précise de ce qu'est l'A.N, et ce pour la simple raison que 
toute tentative de définition utilisera des termes techniques que nous verrons au cours des chapitres 
qui vont suivre. 

Essai de définition: 

L'A.N est l'étude des méthodes constructives permettant d'évaluer numériquement un nombre, une 
fonction ou en général une solution d'un problème donné. 
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6 1. Introduction et définitions 

Par méthode constructives nous entendons les procédures permettant d'obtenir une solution a 
un problème mathématique avec une précision arbitraire en un certain nombre d'étapes. C'est ce 
qu'on appelle Algorithme. 

Le mot "méthodes" qui apparaît dans la définition dépend généralement du problème considéré. 

Définition d'algorithme 

C'est la description précise des opérations successives permettant d'obtenir un certain résultat à 
un problème donné. 

En A.N on s'intéresse à l'étude des conditions d'existence et d'unicité de la solution d'un problème 
ainsi qu'aux performances de l'algorithme telles que la convergence, la stabilité et la précision. 

Exemples. 

a. Soit à multiplier 19 par 17 

19 

17 



Algorithme: < 



133 

19 



323 

b. Soit à effectuer le calcul du P.G.C.D=d de deux entiers relatifs x et y. L'algorithme d'Euclide, 
qui consiste à faire des divisions euclidiennes successives jusqu'à ce que le reste devienne nul donne 
immédiatement le P.G.C.D de x et y. Cet algorithme nous permet aussi de vérifier l'identité de 
Bezout c.à.d qu'il existe u et v de 7L tels que d = xu + yv. 

Algorithme d'Euclide 

x = yq + r Q , < r < y 

y = roqi + ri , < ri < r 

ro = nq 2 + r 2 , < T% < T\ 

ri = r 2 <?3 + ^3 , < T Z < r 2 

r n -2 = r n -iq n -f r„, < r n < r n _i 

r n -l = r n q n+ i + r„+i, < r n+ i < r n 

On arrive forcément à un reste r n+ i = 0, et on conclut que le P.G.C.D de a; et y est le dernier reste 
non nul: P.G.C.D(x, y) = r n . A partir des ces égalités on peut trouver u et v tels que r n = xu + yv 
Programme de calcul du P.G.C.D 

L'algorithme en langage ordinaire s'écrit: 
répéter r := x mod y ; 
x :=y; 
y := r ; 

jusqu'à ce que r = 
(le résultat est le dernier reste non nul) 



program PGCD 
{ calcul du P.G.C.D de deux entiers strictement positifs } 
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1. 2 Les nombres en A.N 7 

var x,y,r: Integer; 

begin writln ('donner x et y (entiers > 0) '); 

repeat Readln(:r,y) Until (x > 0) and (y > 0); 

write('PGCD(',x,7,y,')='); 

repeat r := x mod y 

x :=y; 

y :=r; 

Until r = 0; 

writln(;r) 

End. 

1.2 Les nombres en A.N 

Fréquemment les problèmes d'A.N sont fondés sur des concepts de continûm spatial et temporel. Les 
variables prennent leurs valeurs dans l'ensemble des nombres réels (mécanique, thermodynamique) 
ou des nombres complexes (électricité, électronique, ...). Notez à cet égard qu'un nombre complexe 
du point de vue numérique est la collection de deux nombres réels. 

1.2.1 Connaissance numérique d'un nombre. 

On considérera un réel comme connu numériquement si on connait un certain nombre de chiffres 
de son développement décimal et la précision avec laquelle ce développement approché est donné. 

Exemple. 

y/2= 1.414213 ± HT 6 

7T = 3.14159 ± HT 5 

1.2.2 Connaissance implicite d'un nombre 

Fréquemment la solution d'un problème est un nombre réel ou complexe comme l'illustrent les 
exemples suivants. 

a. Limite d'une suite ou d'une série. 

/— u n + 2 

i. V2 = lim u n tel que u n+ j = — ■ , Vrc G IN, ï/q > 

n-> + OC u n + 1 

ii. 7e = lim (1 H 1 \- . .. -\ log(n)) (je est la constante d'Euler ) 

n^+oo 2 3 n 
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b. Racine d'une équation ou d'un système d'équations. 

• Recherche d'une valeur propre A G (D tel que il existe un X €= IR n : AX = XX et X ^ 0. 

• Recherche d'une solution d'un système: X 6E IR" tel que AX = b A une matrice carrée d'ordre 

n. et 6 E H" 

c. Valeur d'intégrale. 

dt 



it _ r 1 

4"./o ï 



+ t 2 
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d. Optimisation d'une fonction. 

Il s'agit de trouver l'extremum d'une fonction. 

1.3 Les fonctions en A.N 

1.3.1 Connaissance numérique d'une fonction 

Une fonction est connue numériquement si l'on peut calculer "au sens de la connaissance d'un 
nombre réel" sa valeur en tout point de son ensemble de définition. On considérera une fonction 
comme connue numériquement si on peut en calculer une valeur approchée en tout point avec 
indication de la précision de cette valeur. 

Pour les fonctions usuelles "log, sin, cos ... "les algorithmes qui les calculent sont généralement 
fournis par le constructeur de la machine. 

1.3.2 Connaissance implicite d'une fonction 

a. Développement limité au voisinage d'un point. 



i. Log(l - x) = Y^ — 

n=i 



E- 



b. Solution d'équations différentielles. 

• Trouver y tel que y'(x) = f(x,y(x)), pour tout x Ç=]a,b[ avec la condition initiale y(a) = a, a 
donnée. 

• Trouver y tel que y"{x) = f{x,y{x) 1 y'(x)), pour tout x G]a,fe[ avec les conditions au bord 
y(a) = a,y(b) = 0, a et/3 données. 

1.4 Les erreurs en A.N 

Soit x un nombre réel. Soit x une approximation de x. 
On appelle erreur absolue: Ax = x — x. 



1.4.1 Erreurs de méthode. 

Soit à calculer l'intégrale: I = f_ f{t)dt\ f est une fonction continue sur [a, b]. La méthode des 
trapèzes (voir figure) donne pour / la valeur approximative suivante: 

r*=Ar£i(/(a^+i) + /(««)) 

avec k = 2^2 } Xi = a + ik\ i = 0,1, . . .n. 

l'erreur commise en approximant I par T& est: \ T^ — I \. On peut montrer que (voir chapitre 
V): 

| Tfc — I \= o(k ) ordre de l'erreur 

\Tk — I |< Hk . majoration de l'erreur avec H = —(b — a)(max | f"{x) \) x ^[a.b] 
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1. 4 Les erreurs en A.N 9 

1.4.2 Erreurs d'arrondi 

Les opérations de l'algorithme ne sont jamais effectuées exactement. Il y a erreur d'arrondi ou de 
troncature. 

Exemple. 

tt = 0.3333 en simple précision, i = 0.33333333 en double précision. 

1.4.3 Précision machine 

On appelle précision de la machine le réel défini par: 

s = min{e' : é > , 1 + e' > 1} 

avec cette définition on remarque que dans l'ensemble des réels représentés sur une machine 
l'addition n'est pas associative. 

En effet: 

[(1 + s/2) + e/2} machine = 1 alors que 

[1 + (e/2 + e/2)] machine = 1 + e > 1 

1.4.4 Stabilité et perte de précision 

a. Notion de stabilité 

On dira de manière vague qu'un algorithme est numériquement instable si le résultat fi- 
nal est fortement influencé par les erreurs d'arrondi. Le résultat change de la simple à la double 
précision. 

On appelle problème bien conditionné un problème pour lequel on peut trouver un algorithme 
stable pour le résoudre. 

Exemple: 

Le calcul des zéros des polynômes n'est pas toujours bien conditionné. 

b. Perte de précision. 

Nous étudions ce phénomène sur l'exemple suivant: étant donnés deux réels x et y, x et y leurs 
valeurs approchées. 
x = x + Ax , y = y -f- Ay =>■ x + y = x + y + Ax + Ay 

A(x + y) = Ax + Ay 

Si on note © la somme numérique, £ l'erreur relative de méthode commise sur x + y 

x®y = x + y + A(.t + y) = (x + y){l + s) 

l'erreur relative de la méthode numérique est donc donnée par: 

x@y-{x + y) (x + Ax + y + Ay)(l + e) - (x + y) 

x + y x + y 

Ax + A«, , Ax + Ay 

x + y x + y 

où nous avons négligé les erreurs du second ordre dans la dernière équation. 

On peut donc avoir une augmentation significative de l'erreur relative sur la somme de deux réels 

lorsque: 
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10 1. Introduction et définitions 

• (x + y) est très petit devant x et y. 

• — ^ est grand et y -C x 

Exemple. 

Etudions la somme x -\- y + z en supposant que les nombres considérés ont une représentation 
exacte en machine. 
On a x © y = (x 4- y)(l + £i). 

Etudions: a; © y © z: 

(x © y) © z = ((x © y) + z)(l + £ 2 ) = £ + y + ^ + £1(3 + y) + £2(3 + y + 2) + £i£ 2 

L'erreur relative sur cette somme est: S\^—j— + £2- 

On remarque que si x + y + z w (2 ~ — (x + y) ) on a une amplification de l'erreur relative. Dans 

un tel cas, on peut améliorer la précision du résultat en considérant l'un des calculs {x © z) © y ou 

(y © z) © x. 

Numériquement . 

Si x = 1.0000002, y= 1.4999999. 10" 7 et z = -1.0000001 donc x + y + z = 2.4999999.10 -7 . 
Considérons à présent les sommes S± et 5*2 : 

Sj = ( x ®y)®z = (1.0000003) ©2 = 2.10" 7 , S 2 = (x®z)®y = (0.0000001) ©y = 2.4999999. 10" 7 

Les erreurs relatives sur S± et S2 sont définies par: 

_ (x © y) © z - (x + y + z) _ 0.4999999 
651 ~ x + y + z ~~ 2.4999999 ' 

e S2 = 0. 
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